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1  Research  Objectives 

The  classical  asymptotic  description  of  dispersive  pulse  propagation,  initiated  by 
Arnold  Sommerfeld  [1]  and  Leon  Brillouin  [2,  3]  in  1914  using  the  then  newly  de¬ 
veloped  method  of  steepest  descent  due  to  Debye  [4],  and  more  recently  completed 
by  Oughstun  &  Sherman  [5,  6,  7],  Oughstun  [8,  9],  and  Cartwright  &  Oughstun  [10] 
using  modern,  uniform  asymptotic  expansion  techniques,  has  provided  an  accurate 
mathematical  description  of  ultrawideband  pulse  propagation  in  causally  dispersive 
dielectrics  and  conducting  media.  The  accuracy  of  this  asymptotic  description  in¬ 
creases  monotonically  as  the  propagation  distance  z  increases  above  some  character¬ 
istic  distance  Zd  set  by  the  material  dispersion,  typically  given  by  the  e-1  penetration 
depth  at  some  oscillation  frequency  characteristic  of  the  input  pulse.  This  asymp¬ 
totic  description  has  also  resulted  in  a  simple  physical  description  of  dispersive  pulse 
dynamics  [11,  12]  based  upon  the  energy  transport  velocity  [8,  13]  and  attenuation  of 
a  time-harmonic  electromagnetic  plane  wave  in  the  dispersive  medium  that  reduces 
to  the  approximate  group  velocity  description  in  a  specific  limit  of  vanishing  loss  (the 
limit  in  which  the  group  velocity  approximation  is  valid).  What  remains  to  be  done 
in  order  to  complete  this  mathematically  rigorous  theory  of  dispersive  pulse  propa¬ 
gation  is  the  development  of  a  physical  description  of  the  precursor  held  formation 
at  the  molecular  level,  as  it  is  there  that  the  origin  of  the  material  dispersion  occurs 
and  the  precursor  held  formation  first  appears.  Associated  with  this  problem  is  the 
correct  physical  description  of  dispersive  pulse  propagation  in  the  immature  disper¬ 
sion  regime  where  the  precursor  formation  occurs.  Although  it  has  been  asserted  that 
the  group  velocity  description  would  provide  this  near-held  description,  the  research 
conducted  in  this  grant  has  shown  that  this  is  is  not  necessarily  true.  These  two 
related  problems  formed  the  focus  of  this  research  grant. 
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2  Research  Results 


A  molecular  theory  of  optics  was  originally  developed  by  M.  Born  [14],  L.  Rosenfeld 
[15],  and  E.  Lalor  &  E.  Wolf  [16];  a  simpler  asymptotic  description  of  molecular  optics 
may  be  found  in  the  essential  optics  text  [17]  by  M.  Born  &  E.  Wolf.  This  classical 
theory  of  molecular  optics  leads  to  the  fundamentally  important  Ewald-Oseen  extinc¬ 
tion  theorem  which  shows  that  a  time-harmonic  electromagnetic  wave  incident  upon 
a  dielectric  interface  initially  penetrates  into  the  medium  undisturbed,  exciting  the 
molecules  that  comprise  the  medium  which,  in  turn,  produce  a  scattered  wave  held 
that  extinguishes  (through  destructive  interference)  the  incident  plane  wave  held  and 
constructs  the  refracted  plane  wave  held.  This  analysis  thus  reveals  the  fundamental 
physical  processes  involved  in  the  refraction  of  light  at  a  dielectric  interface.  How¬ 
ever,  this  analysis  has  only  been  completed  for  the  special  case  of  a  time-harmonic 
electromagnetic  wave.  The  extension  of  this  theory  of  molecular  optics  to  the  time- 
dependent  case  was  initiated  in  this  grant  with  particular  attention  given  to  the 
ultrawideband  signal,  ultrashort  pulse  regime  so  as  to  reveal  the  fundamental  phys¬ 
ical  processes  involved  in  the  precursor  held  formation  over  the  initial  penetration 
into  the  Lorentz-type  dispersive  half-space.  Because  the  precursor  helds  are  due  to 
the  intricate  interplay  between  the  phase  delay  dispersion  c/vp(co )  =  nr(co)  and  the 
attenuation  dispersion  a(ui)  =  ( u/c)rii(iv )  of  the  dielectric  medium,  where  the  real 
nr(uj)  =  9ft{n(o;)}  and  imaginary  n.i(uj)  =  A{n(u;)}  parts  of  the  complex  index  of 
refraction  n(u)  =  [(e(o;)/eo)  (/ u.//i0 )]^2  are  interrelated  through  the  appropriate  dis¬ 
persion  relations  required  by  causality  [8,  18],  they  are  an  example  of  an  emergent 
physical  process  and,  as  such,  they  present  an  important  example  of  multi-scale  com¬ 
plexity. 


2.1  Integral  Equation  Representation  of  Electromagnetic 
Pulse  Propagation  in  Dispersive  Molecular  Optics 

Consider  the  propagation  of  an  electromagnetic  wave  in  a  homogeneous,  isotropic, 
nonmagnetic  medium  comprised  of  molecules  that  react  to  an  incident  held  like  ideal 
point  dipoles.  The  electric  and  magnetic  held  vectors  E'(r,  t)  and  H'(r,f)  which 
act  on  the  jth  molecular  dipole  in  the  interior  of  the  medium  can  then  be  formally 
separated  into  the  superposition  of  the  incident  electromagnetic  held  vectors  E*(r,  t) 
and  H'(r,f)  that  are  propagating  as  if  they  were  in  vacuum  with  phase  velocity  c 
and  the  contribution  arising  from  all  of  the  other  molecular  dipoles  in  the  medium, 
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so  that 


E'M)  =  E':(r,f)  +  ]TE^(r,f), 

i 

(1) 

H'(r,f)  =  H'(r,t)  +  ^H?,(r,t), 

l 

(2) 

with  j.  At  the  point  ry  where  the  jth  dipole  is  situated,  the  held  of  the  Ith  dipole 

is  given  by  the  pair  of  expressions  [17]  (in  cgs  units) 

Eyw)  =  vkVx^-W, 

Kjt 

(3) 

nyr.c  =  kxM-M, 

(4) 

where  Pt’(t)  is  the  moment  of  the  £th  dipole,  Rjn  =  , 

and  where  the  spatial 

differentiation  is  taken  with  respect  to  the  coordinates  {xv  yv  Zj)  of  the  jth  dipole. 
The  following  simplifying,  but  not  compromising,  approximations  are  now  made: 

1.  first,  the  spatial  distribution  of  the  molecular  dipoles  comprising  the  dielectric 
medium  is  treated  as  continuous  so  that  p i(t)  — >  p(r,f); 

2.  and  second,  the  number  density  Af{ r)  of  dipolar  molecules  is  assumed  to  be  a 
constant  A f  in  the  material  body  and  zero  outside. 


For  a  homogeneous,  isotropic,  locally  linear  dielectric  exhibiting  temporal  dispersion, 
the  total  electric  dipole  moment  per  unit  volume  ( macroscopic  polarization  density ) 
is  given  by  [8] 

P(r,  t)  —  Af  f  a(t  —  t')Fj'(r,t')dt' ,  (5) 

J  —  OO 

where  E'(r,  t)  is  a  spatial  average  of  the  effective  local  electric  held  intensity  E'(r,  t'). 
For  a  causal  medium  response,  Titchmarsh’s  theorem  [8,  18]  requires  that  a(t—t')  =  0 
for  all  t'  >  t,  in  which  case  the  upper  limit  of  integration  may  be  extended  to  +oo. 
The  temporal  Fourier  transform  of  the  resulting  convolution  relation  then  yields 

P(r,  oj)  =  A/”o:(a;)E,(r,a;),  (6) 


where  a(co)  =  f™  a(t)e2UJtdt  is  the  mean  molecular  polarizability  of  the  dielectric 
medium1.  This  molecular  polarizability  characterizes  the  frequency-dependent  linear 
response  of  the  molecules  comprising  the  dielectric  body  to  the  applied  electric  held 

1The  term  “mean”  used  here  indicates  a  spatial  average  over  molecular  sites;  see  §4.4.1  of  [8]. 
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(see  Appendix  D  of  [8]).  Notice  that  the  ‘simple’  linear  relation  expressed  in  Eqs.  (5) 
and  (6)  may  be  generalized  to  include  spatial  inhomogeneity  by  letting  Af  — >  Af( r) 
as  well  as  anisotropy  by  generalizing  the  scalar  quantity  a  to  dyadic  form.  Spatial 
dispersion  may  also  be  included  by  removing  the  condition  imposed  by  spatial  locality 
in  a  generalized  version  of  Eq.  (5)  by  letting  a(t)  — >  d(r  —  r',t  —  t').  For  the  current 
proposal,  only  a  ‘simple’  linear,  causally  dispersive  dielectric  is  considered. 


Upon  going  over  to  a  continuous  distribution  of  the  molecular  dipoles  in  the  dielectric 
medium,  Eqs.  (l)-(2)  become,  with  substitution  from  Eqs.  (3)-(4), 


E'(r,t) 

H'M) 


E*(r,  t)  +  /  V  x  V  x 


a r_ 

R 


a(t  —  t')  E'(r',  t  —  R/c)dt' 


dV, 


H'(r,  t)  +  —  /  V  x 


'Af  f* 

—  /  a(t  —  t  R/c)dt' 

R  . /_™ 


d3r', 


(7) 

(8) 


where  R  =  |r  —  r'|,  and  where  the  indicated  spatial  differentiations  are  taken  with 
respect  to  the  unprimed  coordinates  ( x ,  y,  z )  of  the  held  point  r.  If  the  point  of 
observation  (or  held  point)  r  is  outside  the  dielectric  medium,  then  the  integration  is 
taken  throughout  the  entire  space  occupied  by  the  medium.  However,  if  it  is  inside 
the  medium,  then  a  small  domain  occupied  by  the  molecule  at  that  observation 
point  must  be  excluded  from  the  integration  domain.  As  a  matter  of  course,  one  then 
takes  the  limit  as  T?a  — >  0  after  the  analysis  has  been  completed. 

Equation  (7)  is  an  integro-differential  equation  for  the  electric  held  vector  E'(r,  t). 
When  it  is  solved,  the  magnetic  intensity  vector  H'(r,  t)  associated  with  that  electro¬ 
magnetic  wave  held  is  then  obtained  from  Eq.  (8).  This  pair  of  relations  is  essentially 
equivalent  to  Maxwell’s  equations  for  homogeneous,  isotropic,  locally  linear,  nonmag¬ 
netic  materials  that  exhibit  temporal  dispersion.  Their  solution  has  been  obtained 
[14,  15,  16]  for  the  special  case  when  the  incident  wave  held  is  strictly  monochromatic 
with  angular  frequency  to,  in  which  case  the  mean  polarizability  may  be  treated  as  a 
constant  with  a  =  aS{t  —  t'). 

When  the  incident  electromagnetic  wave  is  pulsed,  and  especially  when  it  is  ultraw- 
ideband  with  respect  to  the  material  dispersion,  special  care  must  be  taken  in  solving 
Eqs.  (7)-(8).  These  integro-differential  equations  are  somewhat  simplified  in  the  tem¬ 
poral  frequency  domain  obtained  by  taking  the  temporal  Fourier- Laplace  transform 
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of  these  relations,  with  the  result 


E'(r,  u)  =  E*(r,u;)  +  Afa(u)  /  V  x  V  x 


E'(r',  u 


ILd 


H'(r, ui)  =  FT(r, u) - Afa{uj)  /  Vx 


E'(r» 


^iuR/c 

1  R 

iujR/c 


d3r', 


R 


dV. 


(9) 


(10) 


Let  E  denote  the  closed  boundary  surface  of  the  dielectric  body  occupying  the  region 
V  G  TZ3.  For  any  observation  point  r  E  T>  inside  the  dielectric,  Eqs.  (6)  and  (9)  may 
be  combined  to  yield  the  relation 


P(r,  u)  =  A fa(u)  E?:(r,  u)  +  Ed(r,  u) 


where 


Ed(r,u;)  =  /  V  x  V  x 


>v-vA 


P(r,o;) 


DiojR/c 

R 


3,„l 


dr 


(11) 


(12) 


denotes  the  contribution  from  the  molecular  dipoles  in  the  dielectric  body  T>  excluding 
the  small  region  about  the  observation  (or  held)  point  r.  The  electric  held  vector 


E'(r,  t)  —  E'(r,  t)  +  Ed(r,  t) 


(13) 


acting  on  the  molecular  dipoles  in  the  interior  of  the  dielectric  body  may  then  be 
obtained  from  the  inverse  Fourier  transform  of  the  expression  [cf.  Eq.  (6)] 

=  (14) 

The  solution  then  depends  upon  the  particular  form  of  the  macroscopic  polarization 
density  P(r,  t).  For  example,  for  the  nondispersive  case  the  dipole  moment  density  is 
taken  to  satisfy  the  homogeneous  wave  equation  V2P(r,  t)  —  (n2/c2) <92P(r,  t)/dt2  =  0 
where  the  constant  n  is  an  unknown  quantity  that  is  to  be  determined  [16,  17]. 

For  a  Lorentz  model  dielectric,  the  macroscopic  polarization  density  P(r,  t)  is  derived 
form  the  microscopic  equation  of  motion 

me  (sj  +  2 5sj  +  WoSj)  =  -qe E' (r,  t),  (15) 

where  me  is  the  mass  of  the  electron  and  qe  the  magnitude  of  the  electronic  charge, 
E((r,f)  being  identified  as  the  effective  local  electric  held  intensity  acting  on  the  jth 
molecular  dipole  as  driving  force  (compare  the  definition  of  the  held  vectors  E'(r,  t) 
and  H'  (r,  t)  in  Eqs.  (1)  and  (2)  here  with  Eq.  (4.200)  of  [8]).  The  additional  force  (in 
cgs  units)  —qe(sj/c)  x  B'(r,  t),  arising  from  the  interaction  of  the  moving  charge  with 
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the  effective  local  magnetic  field,  is  assumed  here  to  be  negligible  in  comparison  to  the 
electric  field  interaction  (due  to  the  smallness  of  the  magnitude  of  this  charge  velocity 
in  comparison  with  the  vacuum  speed  of  light  c),  as  described  in  [23].  Here  cao  is  the 
undamped  resonance  frequency,  and  5  is  the  associated  phenomenological  damping 
constant  of  the  microscopic  oscillator.  The  same  dynamical  equation  of  motion  also 
applies  to  molecular  vibration  modes  when  me  is  replaced  by  the  ionic  mass  and 
is  the  undamped  resonance  frequency  of  the  transverse  vibrational  mode  of  the  ionic 
lattice  structure  [24], 


The  temporal  frequency  transform  of  Eq.  (15)  directly  yields  the  frequency  domain 
solution 


~  /  \  Qe/me  -f-,/  /  \ 

Si(r’w)  =  ^-wg+MuT>’w)' 

(16) 

so  that  the  locally  induced  dipole  moment  Pj(r,  ui)  =  —qeSj(r,uj) 

of  the  jth  molecular 

dipole  is  given  by 

Pj(r,w)  =  “HEj(r,u), 

(17) 

where 

—q2/me 

a  (ca)  =  — ; - ^ - — 

(18) 

to2  —  cOq  +  2  idea 

is  the  mean  molecular  polarizability.  The  macroscopic  polarization  density  is  then 
given  by  the  spatial  average2  of  the  locally  induced  (microscopic)  dipole  moments  as 


P(r,u)  =  M’((pj(r,u)))  =  J\fa(u)((  E'-(r,w))).  (19) 

With  the  identification  that 

E'(r,w)  =  ((E'(r,a;))),  (20) 

the  expression  in  Eq.  (19)  reduces  to  that  given  in  Eq.  (6),  viz. 

P(r,  oj)  =  A/"o:(a;)E,(r,  oj).  (21) 

Following  the  analytical  tack  taken  by  Lalor  and  Wolf  in  [16],  it  is  now  assumed  that 
this  polarization  density  satisfies  the  Helmholtz  equation 

(V2  +  n2(uj)kl)  P(r,  u)  =  0,  (22) 

where  k0  =  co/c  denotes  the  wave  number  in  vacuum,  and  where  the  complex- valued 
quantity  n(ui)  remains  to  be  determined.  As  a  trial  solution,  let  [16] 

P(r,ta)  =  {n2(u)  -  l)  &oQ(r,u;).  (23) 

2The  definition  of  the  spatial  average  of  a  given  quantity  *,  indicated  by  the  double  angle  bracket 
notation  ((*)),  may  be  found  in  §4.1.1  of  Ref.  [8]. 
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The  spectral  quantity  Q(r,  oj)  then  satisfies  the  same  vector  Helmholtz  equation 

(' V2  +  n2(co) fcg)  Q(r,ca)  =  0  (24) 

as  does  the  polarization  density.  The  complex  wave  field  Q(r,  u)  then  “travels”  in  the 
dielectric  body  with  the  complex  phase  velocity  c/n(uj).  In  addition,  it  is  assumed 
that  Q(r,  uj)  has  no  sources3  inside  the  medium,  so  that 


V  •  Q(r,  oj)  =  0 

for  all  r  G  'D.  Combination  of  Eqs.  (11)  and  (23)  then  gives 

1 


Q(r  ,u>)  =  A fa(u>)< 


(n2(u)  —  1)  kt 


2E4(r,cu) 


V  x  V  x 


Iv-V  A 


where 


G(R,u)  = 


0iuR/c 

R 


is  the  free-space  Green’s  function  with  R  =  |r  —  r'|.  With  the  definition  [16] 


Ad(r,w)  =  /  V  x  V  x 
Jv-v  A 


Q(r7,  co)G(R,  ui)  d3r', 


the  relation  in  Eq.  (26)  may  be  expressed  in  a  more  compact  form  as 

1 


Q(r,ta)  =  Afa(co)< 


(n2(w)  -  1  )kl 


E'(r,w)  +  A“(r,w) 


(25) 


Q)r ',u>)G(R,ui)  <fr'y  (26) 


(27) 


(28) 


(29) 


Notice  that  the  mathematical  form  of  the  spectral  field  quantity  Ad(r,u;)  is  precisely 
that  obtained  in  the  classical  description  of  molecular  optics  [16],  which  was  derived 
for  a  monochromatic  wave  field. 


As  the  radius  of  the  small  spherical  region  T> a  surrounding  the  observation  point  r 
shrinks  to  zero,  the  integral  appearing  in  Eq.  (28)  becomes  (see  Appendix  V  of  Born 
and  Wolf  [17]) 


Vx  Vx 


>V-V  A 


Q(r '  ,u)G(R,uj) 


d3r'  =  V  x  V  x 


'V-V  A 


~  87 r  ~ 

Q(r',  uj)G(R,  Ld)d3r' — r-Q(r,  w) 

(30) 


3Note  that  this  macroscopic  polarization  density  is  derived  from  the  spatial  average  of  the  mi¬ 
croscopic  dipoles  induced  by  the  local  effective  electric  field. 
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Because  the  free-space  Green’s  function  G(R ,  u)  describes  a  monochromatic  spherical 
wave  of  angular  frequency  u>  and  wave  number  ko  —  uj/c  in  free-space,  it  then  satisfies 
Helmholtz’  equation 

(V2  +  &q)  G(R,uj)  =  0  (31) 

in  the  region  T>  —  From  Eqs.  (24)  and  (31)  one  finds  that 

QG=(t^IM^v2G-Gv2£!)'  (32) 

so  that,  with  application  of  Green’s  theorem,  the  volume  integral  appearing  on  the 
right-hand  side  of  Eq.  (30)  becomes 


Q(r',  cv)G(R,u)d3r'  = 


'v-vA 


(■ n 2  -  l)fco  Jv-v. A 

1 


QVZG  -  GVZQ)  d3r' 


u"'b  <-"y  l«v 


(n2  —  l)fc(2  )  T*  \  ^  dn'  ~  8 n' 


~dG  <9Q  ,  2  , .  .  . 

alQ9fl-GM|dH.(33) 


where  djdn!  denotes  differentiation  along  the  outward  normal  n!  [not  to  be  confused 
with  the  complex  index  of  refraction  n  =  n(uu)\  to  the  boundary  surface  of  the  region 
T>  —  V a-  Notice  that  a  minus  sign  appears  in  the  second  surface  integral  of  Eq.  (33) 
because  d/dR  is  the  inner  directed  normal  derivative  to  the  spherical  surface  V a- 
The  evaluation  of  the  limit  of  this  inner  surface  integral  over  the  spherical  region 
T> a  as  its  radius  R  approaches  zero  closely  follows  that  given  in  the  derivation  of  the 
integral  theorem  of  Helmholtz  and  Kirchhoff  (see,  for  example,  §8.3.1  of  Born  and 
Wolf  [17])  with  the  result 


£  “  G5s)  dV  = -4^(r)'  (34) 

With  this  substitution,  Eq.  (33)  becomes 

{jf  (C  -  G^)  ^ +4^w}  ■  (®) 

Substitution  of  Eqs.  (30)  and  (35)  into  Eq.  (28)  then  gives 


Ad(r,u;)  = 


77V 


l)kl 


dG  80 

4ttV  x  V  x  Q(r,  oj)  +  V  x  V  x  <f>  I  Q—  -  G-£  ]  dV 


'  dn'  dn' 


(36) 


Because  Q(r,u;)  is  solenoidal  [see  Eq.  (25)],  then,  together  with  Eq.  (24)  it  is  found 
that 

V  x  V  x  Q(r,u;)  =  V  (V  ■  Q(r, u)^  —  V2Q(r, u)  =  n2(w)fcoQ(r,  to),  (37) 
and  Eq.  (36)  becomes 


3  V  n2  -  1 


:VxVx.£  dV- 


Substitution  of  this  expression  into  Eq.  (29)  then  gives 

47t  / n2(cu)  +  2\  .  1  ~ 

1 - T  (  ~TT~\ - 7  )^a  M  Q  r.w 

3  \nz{ui)  —  1/ 


E'M  +  VxV*/  (Qg-G§  dV 


Because  the  left-hand  side  of  this  equation  describes  a  monochromatic  electromagnetic 
wave  with  complex  phase  velocity  c/n{u ),  whereas  the  right-hand  side  describes  a 
wave  that  propagates  with  the  velocity  c,  each  side  must  then  separately  vanish. 
Thus 

(40) 

which  is  known  either  as  the  Lorentz- Lorenz  relation  [25,  26]  or  the  Clausius-Mossotti 
relation  [27,  28],  and 

(41) 

which  is  known  as  the  Ewald-Oseen  extinction  theorem  [29,  30].  Because  the  Lorentz- 
Lorenz  relation  (40)  may  also  be  derived  from  a  separate  line  of  analysis  (see,  for 
example,  §4.4.1  of  [8]),  the  Ewalf-Oseen  extinction  theorem  may  also  be  viewed  as  a 
consequence  of  it  as  applied  to  Eq.  (39). 


The  relation  given  in  Eq.  (41)  expresses  the  extinction  of  the  incident  spectral  wave 
field  component  E*(r,  u)  at  any  point  within  the  dielectric  body  through  destructive 
interference  with  part  of  the  induced  dipole  wave  held.  Notice  that  Eq.  (41)  states 
that  the  extinction  of  this  incident  wave  held  is  brought  about  entirely  by  the  dipoles 
on  the  boundary  surface  of  the  dielectric  body.4  The  extinguished  incident  spectral 

4This  result  is  entirely  consistent  with  the  boundary  conditions  for  the  electric  and  magnetic  field 
vectors  at  an  interface  S  that  are  derived  from  the  macroscopic  Maxwell’s  field  equations. 
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wave  field  at  the  frequency  ca  is  then  replaced  by  the  spectral  wave  field  at  the  same 
frequency  ca  that  is  given  by  [see  Eqs.  (13)  and  (11)] 


E'(r,ca)  =  E*(r,  ca)  +  Ed(r,  ca)  = 


A fa( 


c a 


-P(r,  ca). 


(42) 


With  substitution  from  Eq.  (23)  and  the  Lorentz- Lorenz  relation  (40),  this  frequency- 
domain  expression  for  the  electric  field  vector  acting  on  the  molecular  dipoles  in  the 
dielectric  body  becomes 


E'(r,ca)  = 


n2( ca)  -  l)  fcgQ(r,ca) 


AT  ot{oS) 

^  (n2(u)  +  2)  ^Q(r,w). 


3  (W(ca)  +  2j" 

which  “propagates”  inside  the  dielectric  body  with  the  complex  phase  velocity 

c 


(43) 


vp  = 


n(ca ) 


(44) 


In  terms  of  the  macroscopic  polarization  vector  P(r,  ui),  this  expression  becomes 

4vrn2(ca)  +  2 


E'(r,cu)  = 


3  n2(ta)  —  1 


p(r.w), 


(45) 


where  P(r,ca)  satisfies  the  Helmholtz  equation  (22). 

Finally,  because  D(r, oj)  =  n2(ca)E(r, ta)  in  a  homogeneous,  isotropic,  locally  linear, 
nonmagnetic  medium,  and  D(r,  ta)  =  E(r,ca)  +47rP(r,  ta)  in  a  simple  polarizable  di¬ 
electric,  then  the  elimination  of  the  macroscopic  electric  displacement  vector  spectrum 
from  these  two  expressions  yields 

E(r,ca)  =  P(r,  ca).  (46) 

rr(ca)  —  1 


With  Eq.  (45)  rewritten  as 

E/(r,ca)=  ^  P(r,ca)  +  ^-P(r,ca),  (47) 

rr(ca)  —  1  3 

comparison  with  Eq.  (46)  then  shows  that  the  effective  field  acting  on  the  molecular 
dipoles  in  the  dielectric  body  is  given  by  the  well-known  expression 

E/(r,t)  =  E(r,t)  +  ^P(r,t)  (48) 

after  a  straightforward  Fourier  inversion  into  the  space-time  domain. 
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2.2  Precursor  Field  Formation  in  the  Immature  Dispersion 
Regime 

Based  upon  the  analytical  formulation  of  dispersive  molecular  optics  presented  in  §2.1, 
the  time-domain  development  of  the  precursor  fields  clue  to  an  incident  ultrawideband 
plane  wave  pulse  may  now  be  carefully  investigated.  Associated  with  this  problem 
is  the  correct  physical  description  of  dispersive  pulse  propagation  in  the  immature 
dispersion  regime  where  the  precursor  formation  occurs.  Although  it  is  generally 
asserted  that  the  group  velocity  description  provides  this  near-held  description,  the 
research  conducted  during  this  grant  period  has  shown  that  this  is  not  in  general  true. 


2.2.1  Molecular  Optics  Formulation 

Consider  a  linearly  polarized  plane  wave  pulse  that  is  normally  incident  upon  a  dis¬ 
persive  dielectric  filling  the  positive  half  space  z  >  0.  Let  the  scalar  held  behavior  of 
this  incident  pulse  be  described  by  El( 0~,  t)  =  E0f(t )  where  E0  is  a  constant  and  f(t) 
describes  the  temporal  structure  of  the  pulsed  wave  held.  The  macroscopic  electric 
held  intensity  at  any  penetration  distance  z  >  0  inside  the  dielectric  half-space  in 
then  given  by 

E(z,  oj)  =  2  Eof(ui)e^,  (49) 

1  +  n(u) 

where  the  factor  2/(1  +  n{u))  is  the  Fresnel  transmission  coefficient  [17]  for  normal 
incidence, 

k(oo)  =  —n(u)  (50) 

CO 

is  the  complex  wave  number  in  the  dispersive  dielectric  at  the  frequency  oj,  and  where 

/OO 

f(t)e^dt  (51) 

-CO 

is  the  Fourier- Lap  lace  spectrum  [7,  8]  of  the  initial  pulse  waveform.  The  spectrum  of 
the  induced  macroscopic  polarization  is  then  obtained  from  Eq.  (46)  as 

PM  =  -L(n(a;)  -  l)E0f(oj)Mz  (52) 

for  all  z  >  0,  and  the  spectrum  of  the  electric  held  acting  on  the  molecular  dipoles  in 
the  interior  of  the  dielectric  half-space  is  obtained  from  Eq.  (45)  as 

&(z,u)  =  (53) 

3  n(oj)  +  1 
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for  all  z  >  0.  Finally,  the  temporal  frequency  spectrum  of  the  dipole  held  in  the 
dielectric  half-space  z  >  0  is  obtained  from  Eqs.  (13)  and  (53)  as 

Ed(z,co)  =  \  +  2  EQf(u)eii(u)z  -  E0f(co)eik°z  (54) 

3  n(u>)  +  1 

where  k0  =  oj/c  is  the  wave  number  in  vacuum.  Because  k{uf)  —  k0  =  -(n(ui)  —  1), 
this  expression  can  be  written  in  more  compact  form  as 


E\z,oj)  = 


[2  n2(uj)  +  2  $g(w(a,)_1)z  _ 


3  n(cj)  +  1 


i  Eof^y 


From  Eqs.  (42)  and  (48),  the  propagated  plane  wave  held  in  the  dispersive  dielectric 
occupying  the  half-space  z  >  0  may  be  expressed  as 

E(z,  t )  =  E\z,  t)  +  Ed(z ,  t)  -  jP(Z,  t)  (56) 

This  formulation  will  then  reveal  the  origin  of  the  precursor  helds  through  the  in¬ 
teraction  of  the  dipole  held  Ed(z,t )  with  the  incident  held  El(z,t )  and  the  induced 
polarization  held  P(z,t).  This  comprises  a  major  part  of  my  future  research  focus. 


2.2.2  Energy  and  Group  Velocity  Descriptions 

Both  the  energy  and  group  velocity  descriptions  of  dispersive  pulse  dynamics  are 
based  on  (or  derived  from)  the  exact  Fourier- Laplace  integral  representation  of  plane 
wave  pulse  propagation  [1,  2,  7,  8] 

A(z,t)  =  ±  [  (57) 

where  .4(0,  t)  =  f{t)  describes  the  initial  pulse  at  z  =  0  with  spectrum 

/OO 

me**dt.  (58) 

-OO 

The  complex  phase  function  <j>(uj,9)  appearing  in  Eq.  (57)  is  given  by 

c  r~  i 

4>(u),6)  =  i—  k(u>)z  —  cot  —  iu[n{u)  —  #],  (59) 

where 

kiuj)  =  f3(u)  +  ia(cu)  =  —n(u)  (60) 

c 


12 


is  the  complex  wavenumber  of  the  electromagnetic  wave  in  a  dispersive  nonmagnetic 
medium  characterized  by  the  complex  index  of  refraction  n(u)  =  [ec(o;)/eo]1^2,  and 
where 

6  =  ct/ z  (61) 

is  a  dimensionless  space-time  parameter  defined  for  all  z  >  0.  If  the  complex  index 
of  refraction  is  described  by  a  causal  model,  then  Sommerfeld’s  relativistic  causality 
theorem  [1]  applies,  which  states  that  [7,  9]  if  A(0,t)  —  0  V  t  <  0,  then  A(z,t )  = 
0  V  z  >  0  when  ct/ z  <  1,  in  agreement  with  the  special  theory  of  relativity.  Physically 
meaningful  examples  of  such  causally  dispersive  medium  models  include  the  single 
relaxation  time  r  Debye-model  dielectric 


e(ca)/ e0 


£oo  + 


^sr  ^oo 
1  —  iujT  ’ 


(62) 


with  Coo  >  1  and  esr 
dielectric 


e(0)/eo,  the  single  resonance  frequency  ojq  Lorentz-model 


e(uj)/e0  =  1 


b 2 

co2  —  ujg  +  2  iSu 


(63) 


with  b 2  =  Nq2/me  the  square  of  the  plasma  frequency  with  number  density  N  and 
with  damping  constant  5  >  0,  and  the  Drude  model  conductor 


€c(w)/e0l 


co: 


lo(lo  +  iy) 


(64) 


with  lo2  =  Nq2/me  and  with  damping  constant  7  =  l/rc  given  by  the  inverse  of  the 
relaxation  time  rc  associated  with  the  mean-free  for  free  electrons  in  the  material. 


With  application  of  the  saddle  point  method  of  analysis,  the  mathematical  formula¬ 
tion  given  in  Eqs.  (57)-(60)  leads  to  the  uniform  asymptotic  representations  [9]  (for 
9  >  1) 

A(z,t)  ~  Ab(z,t)  +  Ac(z,t)  (65) 

as  z  — *  00  for  Debye-type  dielectrics, 

A(z,  t )  ~  As(z,  t )  +  Am(z,  t )  +  Ab(z,  t )  +  Ac(z,  t )  (66) 

as  z  — >  00  for  Lorentz-type  dielectrics,  and 

A(z,  t)  ~  As{z,  t)  +  Ab(z,  t )  +  Ac(z,  t)  (67) 

as  z  — >  00  for  Drude- model  conductors.  Here  As(z,  t )  describes  the  first  or  Sommerfeld 
precursor  wave  whose  dynamical  evolution  is  described  by  the  space-time  evolution 
of  the  distant  saddle  points  of  6)  that  appear  in  both  Lorentz-model  dielectrics 
and  Drude- model  conductors,  Ab(z,t )  describes  the  second  or  Brillouin  precursor 
wave  whose  dynamical  evolution  is  described  by  the  space-time  evolution  of  the  near 
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saddle  points  of  cj)(uj,  6)  that  appear  in  both  Debye-  and  Lorentz  model  dielectrics  as 
well  as  in  Drnde-model  conductors,  and  Ac(z,t)  describes  the  signal  contribution  (if 
any).  Uniform  asymptotic  descriptions  of  the  Sommerfeld  precursor  have  been  given 
by  Oughstun  and  Sherman  [6,  7,  9]  for  Lorentz-model  dielectrics  and  by  Cartwright 
and  Oughstun  [31,  9]  for  Drude-model  conductors.  Uniform  asymptotic  descriptions 
of  the  Brillouin  precursor  have  been  given  by  Oughstun  and  Sherman  [6,  7,  9]  and 
by  Cartwright  and  Oughstun  [10]  for  Lorentz-model  dielectrics,  by  Cartwright  and 
Oughstun  [31,  9]  for  Drude-model  conductors,  and  by  Oughstun  [32]  for  Debye-model 
dielectrics.  Finally,  uniform  asymptotic  descriptions  of  the  signal  arrival  and  evolution 
have  been  given  by  Cartwright  and  Oughstun  [10]. 


Figure  1.  Dynamical  evolution  of  a  Heaviside  step-function  signal  in  (left)  a  Lorentz 
model  dielectric  (with  below  resonance  carrier  frequency  u jc  <  u>0)  and  (right)  a 

Debye  model  dielectric. 


A  comparison  of  the  numerically  determined  held  evolution  due  to  a  Heaviside  unit 
step-function  signal 

A(0,  t)  =  fH(t )  =  uH(t)  sin  (uct)  (68) 

where  Unit)  =  0  for  t  <  0  and  Un{t)  =  1  for  f  >  0  with  constant  carrier  frequency  ujc 
is  given  in  Fig.  1  for  Lorentz-  and  Debye-model  dielectrics  at  1,  3,  and  5  absorption 
depths  Zd  [where  Zd  =  a-1  (a ;c)]  shows  the  difference  between  the  two  precursor  held 
structures  in  these  two  dispersive  medium  types.  The  high-frequency  Sommerfeld 
precursor  evolution  is  clearly  visible  immediately  following  the  propagated  wave  front 
(traveling  at  the  vacuum  speed  of  light  c)  in  the  Lorentz-model  dielectric  which  is 
then  followed  by  the  low-frequency  Brillouin  precursor  evolution  and  then  the  main 
signal  evolution  at  the  input  signal  carrier  frequency  u>c,  as  described  by  Eq.  (66) 
with  Am(z,  t)  =  0.5  The  observed  asymmetry  in  both  the  Sommerfeld  and  Brillouin 

5 Because  the  propagation  is  in  a  single  resonance  Lorentz  medium,  the  middle  precursor  field 
Am(z,t)  is  not  present  in  the  propagated  field  evolution. 
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precursor  field  evolutions  is  characteristic  of  the  saddle  point  evolution  in  a  Lorentz 
medium,  as  described  in  [5,  7,  9].  By  comparison,  the  Sommerfeld  precursor  is  com¬ 
pletely  absent  from  the  dynamical  field  evolution  in  a  Debye-model  dielectric.  In 
that  case  the  propagated  field  evolution  begins  with  the  Brillouin  precursor  which 
is  then  followed  by  the  main  signal  evolution  at  the  input  signal  carrier  frequency 
ay,  as  described  by  Eq.  (65).  Notice  that  the  Brillouin  precursor  evolution  is  nearly 
symmetric  about  its  peak  amplitude  point  in  this  case,  a  characteristic  of  the  saddle 
point  evolution  in  a  Debye  medium,  as  described  in  [32,  8].  In  both  the  Lorentz-  and 
Debye-model  cases,  the  peak  amplitude  of  the  Brillouin  precursor  only  decays  alge¬ 
braically  as  z~1/'2  (except  in  the  singular  dispersion  limit  described  below),  this  peak 
amplitude  point  occurring  at  the  space-time  point  6  =  90  =  n(0).  The  accuracy  of  the 
uniform  asymptotic  description  in  describing  each  feature  in  this  type  of  dispersive 
pulse  propagation  phenomena  is  without  peer. 

One  of  the  main  criticism  of  this  asymptotic  theory  is  that  it  relies  upon  absorption 
and  so  cannot  be  applied  to  the  optical  domain  where  the  materials  (optical  quality 
glass,  etc.)  are  weakly  dispersive  and  the  group  velocity  description  supposedly  ap¬ 
plies.  The  material  absorption  in  a  Lorentz  medium  (which  is  entirely  appropriate  for 
the  description  of  the  material  dispersion  in  the  optical  region  of  the  electromagnetic 
spectrum),  described  by  the  attenuation  coefficient  a(u)  =  (u>/c)rii(uj)  for  real  lu, 
decreases  when  either  5  — »  0  or  when  N  — >  0.  This  then  leads  to  the  following  two 
limiting  cases: 


•  Singular  Dispersion  Limit:  In  the  singular  dispersion  limit  the  frequency  dis¬ 
persion  of  n{u)  becomes  increasingly  localized  about  the  resonance  frequency 
ay  as  5  — >  0.  The  term  “singular”  is  used  here  in  the  mathematical  sense  that 
n(u)  fails  to  be  well-behaved  at  a  point  in  some  well-defined  manner,  in  this 
case  in  terms  of  its  differentiability  at  u  =  ay.  All  of  the  material  dispersion 
and  loss  is  then  concentrated  at  the  resonance  frequency. 

•  Weak  Dispersion  Limit:  In  the  weak  dispersion  limit  the  absorption  coefficient 
a  (a;)  vanishes  while  the  material  dispersion  nr(  u)  approaches  unity  at  all  fre¬ 
quencies  to  as  the  molecular  number  density  N  — >  0.  In  this  limit,  dn(co)/du  — »  0 
as  N  — >  0,  in  agreement  with  the  group  velocity  interpretation  of  weak  material 
dispersion. 


The  recent  extension  and  verification  of  the  asymptotic  theory  to  both  the  singular 
and  weak  dispersion  cases  has  shown  that  [22]: 


•  the  peak  amplitude  decay  of  the  Brillouin  precursor  changes  from  z  l-  2  to  z  1/3 
in  the  singular  dispersion  limit  as  6  — >  0,  as  illustrated  in  Fig.  2,  and  that 
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Figure  2.  Numerically  determined  peak  amplitude  decay  (left)  and  slope  of  the  peak 
amplitude  decay  (right)  in  a  Lorentz  model  dielectric  in  the  singular  dispersion  limit 

as  5  — >  0. 


•  a  simple  scaling  relation  between  strong  and  weak  dispersion  cases  is  given  by 


-2 


Ab 


& 


t2  ~  t\  + 


(69) 


the  accuracy  of  this  result  increasing  as  Ni ,  N2  — >  0  with  N2  >  Ni. 


Most  importantly,  these  results  prove  that  the  precursor  fields  persist  in  both  the 
singular  and  weak  dispersion  limits  where  the  group  velocity  description  supposedly 
holds. 


The  energy  velocity  description  [11,  12,  7,  9]  directly  follows  from  the  asymptotic 
description  and  shows  that  the  propagated  wave  field  is  dominated  by  a  single  real 
frequency  Ue  at  each  space-time  point  6  =  ct/z  for  all  z  >  0.  That  frequency  oje  is 
the  angular  frequency  of  the  time-harmonic  wave  with  the  least  attenuation  that  has 
energy  velocity  Ve  equal  to  z/t ,  so  that 


Ve{we)  — 


(70) 


where  [13] 


vE{w)  = 


(71) 


nr(ui)  +  umi(uj)/5 

for  a  single  resonance  Lorentz-model  dielectric.  For  comparison,  the  group  velocity 
in  a  Lorentz-model  dielectric  is  given  by 

,  ,  1  c 


Vg[U)  = 


d/3/doj  nr(u)  +  um'r(oS) 


(72) 
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The  question  then  arises  regarding  what  conditions  are  required  for  ve(w)  — >■  vg(u), 
or  equivalently,  when  does  rii{uj)/8  — >  n'r(uj)7 


For  a  single  resonance  Lorentz-model  dielectric,  the  complex  index  of  refraction 
n{uj)  =  nr{uj)  +  irtiioj )  is  given  by 


with  derivative 


n(  u)  = 
n'(u)  =  b 2 


b2 


\  1/2 
ca2  —  con  +  2/<5ca  / 


cu  +  i8 


\2> 


(73) 


(74) 


n(ca)  (ca2  —  cug  +  2i8uY 

where  b2  =  Nq2/me  is  the  square  of  the  plasma  frequency  with  number  density  N. 
In  the  singular  dispersion  limit  as  5  — >  0, 

b2u 


n'r{oj )  — >■ 


nr( ca)  (ca2  —  ojq) 


2\2  ’ 


(75) 


rii{u) 


->• 


(ca2  —  tag)2  —  62(ta2  —  tag)  +  4<52ca2 
(ta2  —  tag)2  +  4<52ca2 


1/2 


X  (ta2  -  ca02)2  -  b2 (ta2  -  ca2)  +  4<f2ca2  ^  +  ^ 

b2co 

nr(ta)(ta2  —  tag)2  ’ 


(76) 


and  the  energy  and  group  velocities  are  identical  for  all  real  ta.  However,  in  the  weak 
dispersion  limit  as  N  — y  0  (or  equivalently  as  b2  — y  0), 

.  .  (&2/2)(ta2  —  tag)  62ta 

1  J  (ta2  —  tag)2  +  4<52ca2  °  (ta2  —  tag)2  +  4<52ca2  ’ 

v - v- - '  v - v - ' 

nr(u;) 


n'r(  ta)  — )• 


b2oj 

(ta2  —  tag)2  +  452ta2 
^ " 
m(u})/S 


_  1  +  282/{uj2  -  tag) 

l  +  452ta2/(ta2-ta2)2_ 


(78) 


and  the  energy  and  group  velocities  are  approximately  the  same  when  either  ca  <C  cao 
or  ta  3>  tao,  the  approximation  improving  as  8  — >  0. 


This  analysis  then  casts  doubt  on  the  applicability  of  the  group  velocity  method 
to  dispersive  pulse  propagation  phenomena  except  in  the  singular  dispersion  limit. 
A  proper  description  of  dispersive  pulse  propagation  then  needs  to  be  developed  in 
the  immature  dispersion  regime.  Together  with  the  closely  related  problem  of  the 
molecular  theory  of  precursor  field  formation,  this  forms  the  major  focus  of  my  future 
research  direction. 
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3  Peer-Reviewed  Research  Presentations  &  Pub¬ 
lications 


The  following  peer-reviewed  conference  presentations  were  given  during  this  grant 
funding  period: 


1.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “The  Effect  of  Conductivity  on  the 
Brillouin  Precursor,”  2008  USNC/URSI  National  Radio  Science  Meeting ,  Uni¬ 
versity  of  Colorado,  Boulder,  CO  (2008). 

2.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “The  Effect  of  Conductivity  on  the 
Brillouin  Precursor,”  2008  Progress  in  Electromagnetics  Research  Symposium 
(PIERS  2008),  Massachusetts  Institute  of  Technology,  Cambridge,  MA  (2008). 

3.  K.  E.  Oughstun,  “Ultrawideband  Pulse  Propagation  in  Double-Resonance  Lorentz 
Model  Dielectrics,”  2008  Progress  in  Electromagnetics  Research  Symposium 
(PIERS  2008),  Massachusetts  Institute  of  Technology,  Cambridge,  MA  (2008). 

4.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Pulse  Propagation  in  a  Debye  Medium 
with  Static  Conductivity:  The  Search  for  a  Uniform  Expansion,”  Special  Session 
on  Asymptotic  Methods  in  Analysis  with  Applications,  American  Mathematical 
Society  Annual  Meeting,  Washington,  D.C.  (2009). 

5.  K.  E.  Oughstun,  “On  the  Use  &  Application  of  Precursor  Waveforms,”  2009 
13th  International  Symposium  on  Antenna  Technology  and  Applied  Electromag¬ 
netics  and  the  Canadian  Radio  Science  Meeting  (ANTEM/URSI  2009),  Banff 
Conference  Centre,  Banff,  AB,  Canada  (2009). 

6.  K.  E.  Oughstun,  “Dispersive  Pulse  Dynamics  in  the  Few-Cyclc  Pulse  Limit,” 
(Invited  Talk),  Nonlinear  Optics  in  Guided  Geometries,  Weierstrass  Institute 
for  Applied  Analysis  and  Stochastics,  Berlin,  Germany  (2009). 

7.  K.  E.  Oughstun  and  N.  A.  Cartwright,  “Brillouin  Precursor  Decay  in  Dispersive 
Attenuative  Materials,”  2009  IEEE  International  Symposium  on  Antennas  and 
Propagation  and  USNC/URSI  National  Radio  Science  Meeting,  Charleston,  SC 
(2009). 

8.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Pulse  Propagation  in  a  Debye  Medium 
with  Static  Conductivity,”  2009  IEEE  International  Symposium  on  Antennas 
and  Propagation  and  USNC/URSI  National  Radio  Science  Meeting,  Charleston, 
SC  (2009). 
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9.  K.  E.  Oughstun,  “Utility  of  the  Brillouin  Precursor  in  Debye-Type  Dielectrics,” 
2010  USNC/URSI  National  Radio  Science  Meeting,  University  of  Colorado, 
Boulder,  CO  (2010). 

10.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Optical  Precursors  and  the  Group  Ve¬ 
locity  Approximation,”  2010  IEEE  AP-S  International  Symposium  and  USNC/ 
CNC/URSI  National  Radio  Science  Meeting,  Toronto,  Canada  (2010). 

11.  K.  E.  Oughstun,  “Beer’s  Law  and  the  Unique  Penetration  Properties  of  the 
Brillouin  Precursor  in  Complex  Media,”  201 0  IEEE  AP-S  International  Sympo¬ 
sium  and  USNC /CNC/URSI  National  Radio  Science  Meeting,  Toronto,  Canada 
(2010). 


The  following  peer-reviewed  conference  publications  appeared  during  this  grant  fund¬ 
ing  period: 


1.  K.  E.  Oughstun,  “On  the  Use  &  Application  of  Precursor  Waveforms,”  Pro¬ 
ceedings  of  the  2009  13th  International  Symposium  on  Antenna  Technology 
and  Applied  Electromagnetics  and  the  Canadian  Radio  Science  Meeting  (AN- 
TEM/URSI  2009),  pp.  TPlil-4  (2009). 

2.  K.  E.  Oughstun  and  N.  A.  Cartwright,  “Brillouin  Precursor  Decay  in  Dispersive 
Attenuative  Materials,”  Proceedings  of  the  2009  IEEE  International  Symposium 
on  Antennas  and  Propagation  and  USNC/URSI  National  Radio  Science  Meet¬ 
ing,  paper  #528.5  (2009). 

3.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Pulse  Propagation  in  a  Debye  Medium 
with  Static  Conductivity,”  Proceedings  of  the  2009  IEEE  International  Sympo¬ 
sium  on  Antennas  and  Propagation  and  USNC/URSI  National  Radio  Science 
Meeting,  paper  #330.1  (2009). 

4.  K.  E.  Oughstun,  “Beer’s  Law  and  the  Unique  Penetration  Properties  of  the  Bril¬ 
louin  Precursor  in  Complex  Media,”  Proceedings  of  the  2010  IEEE  AP-S  Inter¬ 
national  Symposium  and  USNC/CNC/URSI  National  Radio  Science  Meeting, 
paper  #425.5  (2010). 


The  following  peer-reviewed  journal  publications  appeared  during  this  grant  funding 
period: 


1.  N.  A.  Cartwright  and  K.  E.  Oughstun,  “Ultrawideband  Pulse  Propagation 
through  a  Homogeneous,  Isotropic,  Lossy  Plasma,”  Radio  Science  44,  RS4013- 
RS4024  (2009). 
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2.  C.  L.  Palombini  and  K.  E.  Oughstun,  “Optical  Precursor  Fields  in  Nonlinear 
Pulse  Dynamics,”  Optics  Express  18,  22,  23104-23120  (2010). 

3.  K.  E.  Oughstun,  N.  A.  Cartwright,  D.  J.  Gauthier,  and  H.  Jeong,  “Optical 
Precursors  in  the  Singular  and  Weak  Dispersion  Limits,”  Journal  of  the  Optical 
Society  of  America  B  27,  8,  1664-1670  (2010). 

4.  K.  E.  Oughstun,  N.  A.  Cartwright,  D.  J.  Gauthier,  and  H.  Jeong,  “Optical 
Precursors  in  the  Singular  and  Weak  Dispersion  Limits:  Reply  to  Comment,” 
Journal  of  the  Optical  Society  of  America  B  28,  3,  468-470  (2011). 

In  addition,  the  following  research  monograph  was  also  published: 

K.  E.  Oughstun,  Electromagnetic  &  Optical  Pulse  Propagation.  Volume  2:  Temporal 
Pulse  Dynamics  in  Dispersive,  Attenuative  Media ,  Springer  Series  in  Optical  Sciences 
(Springer,  New  York,  2009). 

Finally,  the  following  seminar  and  workshop  presentations  were  given  during  this 
funding  period: 

1.  “The  Effect  of  Conductivity  on  the  Brillouin  Precursor,”  19th  Annual  AFOSR 
Electromagnetics  Workshop,  San  Antonio,  TX  (January  8-10,  2008). 

2.  “On  the  Use  &  Application  of  Precursor  Waveforms,”  20th  Annual  AFOSR 
Electromagnetics  Workshop,  San  Antonio,  TX  (January  6-8,  2009). 

3.  “Dispersive  Pulse  Propagation:  History,  Theory,  &  Applications,”  Engineering 
Seminar,  University  of  Toronto  (February  12,  2009). 

4.  “Dispersive  Pulse  Propagation:  History,  Theory,  &  Applications,”  Department 
of  Physics  Seminar,  Naval  Postgraduate  School,  Monterey,  CA  (March  27, 
2009). 

5.  “Precursor  Field  Behavior:  Final  Thoughts,”  21st  Annual  AFOSR  Electromag¬ 
netics  Workshop,  San  Antonio,  TX  (January  5-7,  2010). 

6.  “Precursor  Wave  Fields:  The  Characteristic  Wave  Structure  of  Dispersive  At¬ 
tenuative  Media,”  Institute  of  Optics  Colloquium,  The  Institute  of  Optics,  The 
University  of  Rochester  (April  12,  2010). 

7.  “Multus  Tumultus  de  Nihilo,”  22nd  Annual  AFOSR  Electromagnetics  Work¬ 
shop,  San  Antonio,  TX  (January  4-6,  2010). 

The  Latin  phrase  “Multus  Tumultus  de  Nihilo”  used  as  the  title  in  this  last  workshop 
presentation  was  a  reflection  on  the  “Reply  to  Comment”  paper  listed  above  that 
formed  the  focus  of  this  presentation. 
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